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Abstract. We develop a simple three compartment model based on mass 
balance equations which quantitatively describes the dynamics of breath methane 
concentration profiles during exercise on an ergometer. With the help of this model it 
is possible to estimate the endogenous production rate of methane in the large intestine 
by measuring breath gas concentrations of methane. 
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Various studies demonstrated that methane in humans majorly originates from 
anaerobic fermentation by methanogens in the large intestine. Methane can then 
traverse the intestinal mucosa and be absorbed into the systemic circulation. Since 
methane has a low solubility in blood, it is rapidly excreted by the lungs. It is a 
generally accepted criterion that a subject is considered to be a methane producer if the 
methane concentration in exhaled breath exceeds the ambient air level by 1 ppm EE]- 
Approximately 30-50% of adults were found to be methane producers [30] . Considering 
methane production, gender, age, and ethnic differences were observed [2111211 EH 125] . 
Additionally, a signihcant day-to-day variation was reported I2D1- However, the factors 
influencing the number of methanogens and the amount of methane produced are still 
unexplored. 

The interaction between methanogens and gut function is an extensively studied 
held. Breath methane tests and culture based methods have traditionally been used to 
characterize methanogen populations [7]. Culture based methods have high sensitivity; 
however they are cumbersome and time-consuming. Nevertheless the methane breath 
test is a convenient, quick and effective method for the assessment of methanogen 
populations; therefore it is increasingly used in the diagnostics of certain gastrointestinal 
conditions. In clinical practice, a combined hydrogen and methane breath test has been 
shown to be superior for the diagnosis of carbohydrate malabsorption syndromes and 
small intestinal bacterial overgrowth [7] . It is commonly accepted that breath methane is 
associated with alterations in intestinal motility, and it is strictly related to constipation 
[23 EH El]. Additionally, numerous studies have found correlations between breath 
methane levels and diseases including colon-rectal cancer, irritable bowel syndrome and 
inflammatory bowel disease [27] [281 EH EH- However, the results are controversial and 
the impact of endogenous bacterial methane generation on health is still not known with 
certainty. 

Although numerous studies have conducted methane breath tests, there are only 
relatively few studies that investigated the routes of methane excretion, i.e., the 
correlation between methane concentration in breath and in the gut [3]. It is generally 
assumed that methane is not utilized by humans, and approximately 20% of the methane 
produced by anaerobic fermentation is excreted by breath. The remaining 80% is lost 
by flatus [3]. 

It is worthwhile to note that a recent paper by Boros et ah reviewed the possible 
role of methane as a gasotransmitter [3] . It provides some evidence with respect to non- 
bacterial generation of methane in target cells which is possibly linked to mitochondrial 
dysfunction. Furthermore, methane-rich saline is hypothesized of having an anti- 
oxidative effect [5]. 

Breath tests can be performed even in real time allowing to monitor biological 
processes in the body. In our recent study the dynamics of endogenous methane release 
through the respiratory system have been investigated by measuring breath methane 
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concentration profiles during exercise on an ergometer [29]. The qualitative behavior 
of such prohles was in good agreement with the Farhi equation |9] but the quantitative 
behavior deviated. The aim of this article is to develop a simple three compartment 
model to describe and explain quantitatively the observed breath methane concentration 
prohles. The present model can serve as a tool to estimate the endogenous production 
rate of methane in the large intestine from exhaled methane concentrations. 

A list of symbols used is provided in Appendix A. 

2. Measurements 
2.1. Setup 

End-tidal methane concentration prohles were obtained by means of a real-time setup 
designed for synchronized measurements of exhaled breath VOCs as well as a number of 
respiratory and hemodynamic parameters. Our instrumentation has successfully been 
applied for gathering continuous data streams of these quantities during ergometer 
challenges [H] as well as during sleep studies [13]. These investigations aimed at 
evaluating the impact of breathing patterns, cardiac output or blood pressure on the 
observed breath concentration and permit a thorough study of characteristic changes in 
VOC output following variations in ventilation or perfusion. An extensive description 
of the technical details is given in a previous work [T3]. 

In brief, the core of the mentioned setup consists of a head mask spirometer 
system allowing for the standardized extraction of arbitrary exhalation segments, 
which subsequently are directed into a Selective Reagent Ionization Proton Transfer 
Reaction Time of Flight Mass Spectrometer (SRI-PTR-TOF-MS, lonicon Analytik 
GmbH, Innsbruck, Austria) for online analysis. This analytical technique has proven to 
be a sensitive method for the quantification of volatile molecular species M down to the 
ppb (parts per billion) range. To measure methane we took advantage of the reaction 
of the primary precursor with methane [U El [ 33 ] 

+ C H 4 ^ C Ha OOH+ + H. 

Count rates of the resulting product ion appear at the specihed mass-to-charge ratio 
mfz = 47.0128 (see Figure and hgure 6 in [13]) and can subsequently be converted 
to absolute concentrations by means of calibrations factors obtained from analyzing 
calibrations mixtures containing a known amount of methane and humidity. 

So far, only some preliminary measurements were carried out by means of the 
setup described above. Two healthy methane producing adult volunteers (one male, 
one female) were asked to perform several ergometer challenges of approximately 6 
minutes rest, 17 minutes with 75 Watts, and then approximately 6 minutes rest again. 
The exact protocol was: 

• seconds 0-380: the volunteer rests on the ergometer 

• seconds 380-1400: the volunteer pedals at a constant workload of 75 Watts 
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Figure 1. Spectrum of methane as measured by SRI-PTR-TOF-MS using Oj primary 
ions. 


• seconds 1400-1800: the volunteer rests on the ergometer 

Figure shows a tyical result of such an ergometer session for one volunteer. While the 
number of probands is certainly very limited, the relative changes of breath methane 
concentrations are in good agreement with similar measurements employing a different 
analytical set up as described in a recent work [29] (see figure 1 therein). 


3. Modeling methane distribution in the body 

3.1. Methane exchange in the lungs 

In humans, methane is mainly produced by enteric bacteria in the large intestine and 
distributed within the body by the venous blood leaving the intestine. When reaching 
the lungs, it is partially released into breath. The amount of methane transported at 
time t to and from the lungs via blood flow is given by 

Qcmcfit)-cfit)), 

where Qc denotes the cardiac output, the averaged mixed venous concentration, and 
Ca is the arterial concentration. 

On the other hand one in- and exhales the amount 

I4(t)(C'i-C'A(t)), 

where Va denotes the alveolar ventilation, Ci denotes the concentration in inhaled air, 
and Ca the alveolar concentration. While Ci is assumed to be zero for many endogenous 
VOCs, the current average atmospheric methane concentration level is about 1.8 [ppm] 
and can hence not be neglectec^ 


f Typical room air concentrations are often even higher than 1.8 ppm. 
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Figure 2. Typical result of an ergometer session for one single volunteer. Average 
values: cardiac output (green), rest: Qc = 5.41 [£/min]; 75 Watts: Qc = 11.07 [£/niin]; 
alveolar ventilation (red), rest: Va = 10.69 [f/min]; 75 Watts: Va = 33.12 [.^/min]; 
and exhaled end-tidal (nose sampling) methane levels (blue), rest: Ca = 31.08 [ppm]; 
75 Watts: Ca = 11.92 [ppm], room air concentration of methane: 3.37 [ppm]. 


Combining these two terms leads to the following mass balance equation for the 
lung^ 

Va^ = Va{Ci - Ca) + 4(a - a), (1) 

where Va denotes the volume of the lung. Both sides of Equation ([^ have units 
yumol/min (compare Appendix A). 

If the system is in an equilibrium state (e.g., stationary at rest) Equation ([^ reads 
0 = Va{Ci — C'a(C'i)) + QfiCfiCi) — Cfi) and using Henry’s law 

Ca = Ab:air Ca (2) 


we obtain 


Ca(C'i) 


Ci ^ CfiCi) 

+ 1 •^b:air + T 


(3) 


where r = Va/Qc is the ventilation-perfusion ratio and Ab:air denotes the blood:air 
partition coefficient. 

Remark: The modeling approach followed above is only valid for VOCs with 
blood:air partition coefficient less than 10, i.e., compounds for which the upper airways 


§ For notational convenience we have dropped the time variable t, i.e., we write Cx instead of Cx{t), 
etc. Cx denotes the instant or averaged concentration of X over a small sampling period r, i.e., 
Cx{t) = 1/t Cx{s)ds. 


Qc, Vk [1/: 
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have no influence on the observable breath concentrations j2]. Methane with a blood:air 
partition coefficient Abiair = 0.066 [26] fulfills this requirement. 

Since Ab:air for methane is so small we get Cv(Ci) ~ C'v(O), 1 + ~ 1, and 

Ab;air + r ^ r. From this follows that for methane it suffices to subtract the inhaled 
methane concentration to correct for room air levels (see a previous work for more 
details [31]). 

Thus Equation (|^ can be simplified to 

c'a(o) = c'a(c'i) - a = 2^ a(o) = - a(o). (4) 

14 r 

When a subject is under constant conditions at rest, Cy is approximately constant. 
From Equation (|^ it may then be concluded that variations in the alveolar concentration 
Ca(0) = (C'a(Ci) — Cl) directly reflect changes in ventilation (e.g., due to altered 
breathing frequency) and perfusion (e.g., due to altered heart rate). This can be tested 
by forced hypo- and hyperventilation at rest as shown in figure 5 in a previous work 

ra. 

The ventilation-perfusion ratio r is approximately one at rest but substantially 
increases for a moderate exercise regime at 75 Watts, since the cardiac output 
increases approximately two-fold while the ventilation increases three- to four-fold [H] . 
Consequently, one would expect from Equation (|^ that the alveolar methane 
concentration should decrease by a factor of approximately 1.5-2 when exercising at 
that workload. 

Contrary to this prediction, measurements of breath methane concentrations show 
a drop by a factor of 3 to 4 when exercising at 75 Watts [151121], see also Figure]^ 

3.2. A three compartment model 

The intuitive rationale for this phenomenon is as follows. The intestinal bacteria are the 
main source of methane. At rest, the intestine receives about 15% of the total blood flow 
of approximately 5 £/min, leading to an absolute perfusion of approximately 0.75 £/min, 
which is matched to the metabolic needs of gut tissue. When exercising moderately, this 
absolute blood flow to the intestine may be assumed constant, since its metabolic needs 
remain largely unchanged. However, the relative (fractional) blood flow to the intestine 
decreases, as a major part of total cardiac output is now directed to the working muscles. 
As a result, the relative contribution of intestinal venous blood (characterized by high 
methane concentrations) to mixed venous blood will be reduced, causing the mixed 
venous methane concentration to drop. The decrease in breath methane concentrations 
during exercise may hence be interpreted as a combination of two separate effects: 
an increased dilution within the lungs due to an increased ventilation-perfusion-ratio 
(cf. Equation (|^) and an additional reduction of the mixed venous concentration levels 
due to a reduced fractional perfusion of the intestine. 

In order to mathematically capture the mechanism illustrated above, we developed a 
three compartment model based on mass balance equations, similar to previous modeling 
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efforts, e.g., with respect to isoprene The model consists of a lung compartment, a 
gut compartment (intestine), and a richly perfused compartment which comprises the 
rest of the body as shown in Figure 


a 


a 




gut, 15*" 


Va 


a 


dJ 


Qc 


a. 


(l—ggutlQc (~i 
O'rptjb^ C/a 




rpt 


Li-pt 

""met 


^gutQc 


a 






lung compartment 


richly perfused 
compartment 


gut compartment 


Figure 3. Three compartment model for methane: lung compartment with gas 
exchange, gut compartment with production of methane by enteric bacteria, and 
richly perfused tissue compartment containing the rest of the body including muscles 
(possible but small production and metabolic rate) 


The mass balance equation for the lung compartment has already been derived in 
Equation ([^. Arterial blood leaving the heart with concentration Ca. is divided into two 
blood flows O'gutQc and (1 — qgnt)Qc, where ggut denotes the fractional blood flow to the 
intestine. 

The molar flow to and from the gut compartment is given by qgntQcCa and 
5'gutQcAb:gutCgut, respectively, where the proportional factor Abigut is the corresponding 
blood:tissue partition coefficient. This yields the following mass balance equation for 
the gut compartment (intestine): 
dC 

= q,utQc{Ca - Ab^gutC'gut) + h k^\ (5) 

Here, Hgut denotes the effective volum^of the gut. The factor /i 0.2 respects the fact 
that 80% of methane is lost by flatus and therefore does not enter the blood stream [3]. 
In addition we assume that within the time frame of the ergometer sessions presented, 
the net production rate of methane stays constant and a possible metabolization in 
the large intestine can be respected by a correction of Both sides of Equation (|^ 
have units /xmol/min (compare Appendix A). 

II The vascular blood compartment and the intracellular tissue compartment are assumed to be in an 
equilibrium and therefore can be combined into one single gut compartment with an effective volume. 
For more details about effective volumes compare appendix 2 in a previous paper m- 
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Analogously, for the richly perfused tissue compartment containing the rest of the 
body including muscles we get 
dC 

rpt _ n „ .r. Ab^pt^pt) - CetAb:rptapt + , 




— (1 ~ 5'gut)Qc(Ca 


“pr 


( 6 ) 


dt 

where V^pt denotes the effective volume of this compartment, respects a possible 
small nonbacterial production rate and represents a possible small metabolic rat^ 
Both sides of Equation ([^ have units /imol/min (compare Appendix A). 

Remark: According to Bond [3] both and k'^^ are very small and hence can be 
neglected in a hrst modeling approach. 

The mixed venous concentration is given by the weighted sum of the two body 
compartment concentrations 

C*v (1 Q'gut)Ab:rptC*rpt T Q'gutAbigutC'gut- (7) 

The total mass balance given by the Equations ([^, (|^, and ([^ constitutes a coupled 
system of three hrst order ordinary differential equations (ODEs) of the form 


^c(f) = g(f, c{t)) =: A{f) c{t) + b(f) 
for the three unknown concentrations 

c{t) = {cfit),a,fit),c,^fit)). 

The matrix A{t) and the vector h{t) are given by 


( 8 ) 


(9) 


A{t) = 


Va 

(1-fe,.(«))» 

V 9 sn(*)^ 


Va 

TP 

met 


ggut(^)Ab:airAb:gut%^ \ 


(l-ggut(t))Ab:rptQc(t) + Ab:rptfcme1 

Vrpt 
0 


Va 


0 


Agutit) Ab:gut 


/ \ VaC*)^ \ 


b(f) = 


V 


Va 

l.rpt 

/vpr 

pgut 

^Vgut 


( 10 ) 


/ 


All external inputs (yA(t),Qc(t),Ci) affecting the system can be measured by 
means of the experimental setup and are therefore assumed to be known. The 
partition coefficients Ab:air, Ab:rpt, Abigut may partially be derived from literature values, 
see Section [3^ 

Model parameters that are a priori unknown and not directly measurable include the 
metabolic rate fh® production rates and k^^, as well as the effective volumes 
f^pt, V, which inhuence the time constants for achieving a steady state. These 
will either have to be hxed at best-guess values or estimated from the measurement data 


by means of a suitable parameter estimation scheme, see Section 3.4 


^ Here we used the usual convention to multiply kfffi by Ab:rpt- It would be more natural to use 
only. 
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As explained in the model rationale, the absolnte blood flow throngh the intestine 
is postnlated to stay approximately constant dnring moderate exercise. We therefore 
nse the following simple model for the fractional blood flow 

qgnt{t) = go = 0.15 (11) 

where Qc,rest is the average total blood flow (cardiac ontpnt) at rest. 

In addition, the methane concentration in exhaled end-tidal air is measnred and 
assnmed to be eqnal to the alveolar concentration 

y{t) := C^eUt) = C'a(^) = Ab;iirC'a(^)- (12) 


3.3. Steady state analysis 


When in a steady state the system of differential eqnations rednces to the following 
simple linear algebraic system 

0 = Va{Ci-Ca) + Qc{C,-C^), 

0 qgutQci^Ca '^bigutC'gut) A /^^pr ; 

0 = (1 — Q'gut)Qc(C*a ~ Ab:rptC*rpt) ~ ^riLt-^birptCrpt + ^p?*- (13) 

Solving with respect to CgutjC'rpt, and yields 


Crnt — 


Ab:air {I - q^ut)Qc CA + 


Ab:rpt (1 — q'gut)Qc + ^met 

Ab:air ~ Cj) (1 — Ab:air(l ~ ?gut)Qc Ca + k^ffi 


Cgnt - . 

Ah: 


b:gut Qgnt 

KT =-(MCa-Cj) + 


1 ^ ^ N (1 qgnt)Q^ 


/i 


(1 ~ ?gut)Qc + ^met 


Ab:gut Qgut (1 — qgut)Qc + ^mel 

(Ab:airCe\C^A-fc;f 


(14) 


If we assnme the nonbacterial prodnction and the metabolic rate in the richly 
perfnsed compartment to be negligible we set = 0 and = 0, respectively. Then 
Eqnation (14) simplifies to 

C'rpt(C'l) =^C'a(C'i), 

'^birpt 

c'gat(a) = ^ c'a(q) + —(c'A(a) - a), 




Ab:gut 
1 .. 


5'gut Ab: 


gut 


pr 


= -VAiCAiCi)-Ci). 
S' 


( 15 ) 


Fnrthermore, we recall that 


C*v(Cl) — (1 ~ <?gut)Ab;air Ca(Ci) + <?gutAb:gut C*gut(Cl), 

Ca(Ci) = Ab:aii.CA(Ci). 
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Here we explicitely indicated the dependence of the various quantities on the inhaled 
concentration Cj. 

From Equation (15) we conclude that for a steady state (e.g., at rest or at a 
moderate constant workload): 

(i) The methane concentration Crpt(C'i) in the richly perfused tissue compartment is 
proportional to the alveolar concentration. However, C'rpt(C'i) is much smaller than 
C'a(Ci) since Ab:air is very small. 

(ii) Analogously, since Ab:air is small for Cgut we get 


C,ut{Ci) 
or, vice versa, 
C'a(O) ^ 


Qgnt Ab:gut 


(C'a(C'i) - Q) = 


Qgut Ab;gut 


C'a(O) ^ C'g.tW 


Q'gut Ab:gut 


C'gut(O), 


(16) 


111 


showing that the breath methane concentration is roughly proportional to the 
fractional intestinal blood flow ^gut- 

Since we expect to be constant on a “medium time scale” (e.g., during an 
ergometer session) we obtain 


CA{0)=CA{Cj)-Ci = fikl-^^ 




(17) 


Thus the product C'a(O) x Va does not change when switching from one stationary 
regime to another, e.g. when switching from a resting steady state to an exercise 
steady state at 75 W, viz.. 


— VA,restCA,rest{0) — /c®”* — — VAJdWatts CAJ5Watts(0) ■ 

/i ^ h 


(18) 


This explains the experimental findings of a recent work [21] (see figure 3 therein). 
The production rate of methane in the intestine can therefore be estimated by 
taking the product of average steady state values of Va and Ca(0), 


k^f = -VACAiO). 

/i 


(19) 


3 . 4 . Simulation of an ergometer session and parameter estimation 

In this section we calibrate the proposed model based on the physiological data presented 
in Figure corresponding to one single representative volunteer following the line of 
a previous work [12]. It will turn out that the model appears to be flexible enough 
to capture the methane profiles in exhaled breath generally observed during moderate 
workload ergometer challenges as conducted in a recent work [29]. In a first attempt 
we set the parameter describing a possible small nonbacterial production rate to zero, 
i.e., we fix fcpP* = 0. For VVpt and Va we use the nominal values Va = 4.1 [4] and 
Vj.pt = 15.22 [£] (compare with table C.l 1 in a previous work [12]), and V^ut = 1 [^]- 
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The remaining unspecified parameters pj G be estimated from 

the knowledge of measured breath methane concentrations y by means of parameter 
estimation. More specifically, the subject-dependent parameter vector 


P = ^met) 

as well as the nominal steady state levels Cq = c{to) can be extracted by solving the 
ordinary least squares problem 


argmin (pi - CA{ti)y 

P, CO ^^ 


( 20 ) 


i=0 


( 21 ) 

, and Pi = Cmeas.i is the measured 


subject to the constraints 

g(fo, Co, p) = 0 (steady state) 

P, Co > 0 (positivity). 

Here, g is the right-hand side of the ODE system 
end-tidal methane concentration at time instant ti {to = 0). 

For this purpose the measured physiological functions Va and Qc were converted 
to input function handles by applying a local smoothing procedure to the associated 
data and interpolating the resulting profiles with splines. Furthermore, while the richly 
perfused compartment so far has been treated as an abstract control volume without 
particular reference to any specific tissue group, for identifiability reasons we now set 
Ab:rpt = 1 as well as Ab:gut = 1 which corresponds to the in vitro blooddissue methane 
partition coefficient for brain tissue in rabbits [22], as currently no further values have 
been published. Initial concentrations and fitted parameters are given in Table 


Variable 

Symbol 

Fitted value (units) 

Production intestine 

kp 

51.4 [pmol/min] 

Metabolic rate 

iLrpt 

^met 

0.01 [^/min] 

Initial concentration alveoli {t = 0) 

Ca 

1.15 [yvaol/t] 

Initial concentration rpt (t = 0) 

C'rpt 

0.076 [yuYo\/£] 

Initial concentration intestine (t = 0) 

Cgut 

13.6 [^mol/£] 


Table 1. Decisive model parameters resulting from the fit in Fig. 


All estimated quantities for the test subject under scrutiny take values in a 
physiologically plausible range. According to Equations (|^ and ([^, arterial and mixed 
venous blood concentrations at the start of the experiment are estimated for f = 0 as 
Ca = 0.076 fimo\/£ and = 2.1 jj,mol/£, respectively. Total endogenous production 
is estimated to equal approximately 51.4 [/xmol/min]. The simulation indicates also a 
very small metabolic rate of 0.01 [£/min], which is negligible compared to the production 
rate. 

The results of the simulation are presented in Fig. The first panel of Figure 
shows that the methane concentration profiles obtained from the experiment and from 
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the model are in good agreement. This suggests that the three-compartment model 
can describe quantitatively the methane prohle changes during an ergometer challenge, 
while the Farhi equation provided solely qualitative agreement [2U] . 

4. Conclusion 

Despite the fact that methane breath tests are now widely accepted in clinical practice, 
a quantitative description of the routes of methane excretion is still lacking. The 
present paper intends to £11 this gap by introducing a model for the distribution of 
methane in various parts of the human body. Particularly, we aimed at capturing the 
exhalation kinetics of breath methane in response to exercise. Classical pulmonary inert 
gas elimination theory according to the Farhi equation [HI is deficient in this context, 
as the experimentally observed drop of breath methane concentrations during moderate 
exercise cannot be explained by an altered pulmonary excretion alone. Apart from an 
increased dilution of breath methane within the lungs (due to a rise in the ventilation- 
perfusion ratio r), exercise will also alter the fractional (but not the absolute) perfusion 
of the intestine, which represents the major production site of methane in the body. 
This in turn leads to an additional reduction of mixed venous methane concentrations. 
On the basis of this rationale, a three compartment model extending the original 
Farhi formalism was developed and demonstrated to be in excellent agreement with 
measurement data obtained from a previous study as well as from a SRI-PTR-TOF-MS 
setup presented in this paper. 

From the model equations it can be deduced that under constant resting or 
workload conditions the breath methane concentration Oa(0) is affected by changes 
of the ventilation-perfusion ratio r but also by changes of the fractional intestinal blood 
flow ggut, viz., 

C'a(O) ~ ^Ab,gutC'g„t(0). (22) 

This equation provides a mechanistic physiological rationale for explaining a part of the 
substantial intra-subject variability commonly observed in methane breath tests |H^ITU] . 
In particular, alveolar ventilation can change considerably during breath sampling, since 
patients tend to hyperventilate in such a situation |6|. In this context, it has been 
suggested to normalize breath methane concentrations with respect to CO 2 levels in 
order to improve the repeatability of breath measurements from the same individual [19] . 
Alternatively, as follows from 

CA{0)VA = fik^T, (23) 

the present model points towards l -®-5 fl^® inverse of alveolar ventilation, as an 

appropriate normalization factor for steady-state breath methane concentrations, as this 
allows for a direct estimation of the underlying endogenous methane production rate 
in the intestine. Here, pi ~ 0.2 is a constant factor reflecting the expected methane 
loss due to flatus. For perspective, taking the average resting values from Figure 



Modeling of breath methane concentration profiles 


13 






Figure 4. First panel: simulation of end-tidal methane concentration behavior 
during exercise conditions, cf. Fig. Second panel: predicted methane concentrations 
in mixed venous blood (Cv). Third panel: venous blood concentration returning from 
the gut (Cgut) and returning from the richly perfused tissue (Crpt). Fourth panel: 
predicted profile of the fractional gut blood flow qg^t according to Equation (11). 
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Va = 10.69 [£/min] and C'a(O) = (31.08 —3.37)/27 [p,mo\/i] yields an estimated (resting) 
prodnction rate of 54.9 [/rmol/min]. Analogously, taking average values for a workload 
of 75 W, V A = 33.12 [£/min] and C'a(O) = (11.92 —3.37)/27 [/imol/f'] yields an estimated 
(workload) production rate of 52.4 [/imol/min]. Both estimates are in good agreement 
with the value obtained from htting the model dynamics to the data, see Table [1} In 
particular, note that the estimated endogenous methane production rate during rest and 
exercise is roughly constant (which is in accordance with physiological intuition), while 
the average breath methane concentrations during these two phases differ by a factor 
of roughly 2.6. This proves the efficiency of the above-mentioned normalization scheme 
with respect to reducing the inherent physiological variability due to Equation (22). 

In this sense, the model is expected to contribute towards an improved 
comparability between breath methane measurements as well as towards a better 
quantitative understanding of the correlation between exhaled methane and gut methane 
production in general. Measuring breath methane in combination with the present three 
compartment model can serve as a useful tool to assess endogenous methane production, 
the latter being associated with several gastrointestinal dysfunctions. 
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Parameter 

Symbol 

Unit 

cardiac output 

Qc 


alveolar ventilation 

Va 

[£/niin] 

ventilation-perfusion ratio 

r 

[1] 

averaged mixed venous concentration 

a 

[/imol/£] 

arterial concentration 

a 

[pmol/£] 

inhaled air concentration 

a 

[ppm] 

alveolar air concentration 

C'a 

[ppm] 

richly perfused compartment concentration 

C*rpt 

[pmol/£] 

gut compartment concentration 

Cgut 

[pmol/£] 

measured exhaled concentration 

^meas 

[ppm] 

lung volume 

14 


effective volume of the richly perfused compartment 

Kpt 


effective volume of the gut compartment 

Vgut 


metabolic rate in the richly perfused compartment 

Li-pt 

""met 

[f'/iiiin] 

production rate in the richly perfused compartment 


[pmol / min] 

production rate in the gut compartment 


[pmol/min] 

blood:air partition coefficient 

'^b:air 

[1] 

blood:richly perfused compartment partition coefficient 

-^b:rpt 

[1] 

blood:gut compartment partition coefficient 

-^b:gut 

[1] 

fractional blood flow to the intestine 

?gut 

[1] 

fractional loss of methane due to flatus 

fi 

[1] 


Conversion from [ppb] to [nmol/i]: 


A concentration of x [ppb] corresponds to [nmol/t'j (alternatively, x [ppm] correspond 
to ^ [pmol/t']). The molar volume Idn can be derived from the ideal gas equation (which 
can safely be used for trace gases). For a measured pressure p of 94600 [Pa] and a breath 
temperature of 34 [°C] we get 


Frrt 


R T 

p 


8.314472 (273.15 + 34) 
94600 


27 [£'. 
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